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ABSTRACT 

The  diffraction  of  a  weak  hydromagnetic  shock  wave  by 
a  rigid  perfectly  conducting  half  plane  is  considered.   The 
medi^um  is  assixmed  to  be  a  perfect  compressible  fluid  with 
perfect  electrical  conductivity,  permeated  by  a  \iniforra 
magnetic  field.   The  medium  occupies  the  space  outside  a 
rigid  perfectly  conducting  half  plane.   The  reflected  and 
diffracted  waves  resulting  after  an  incident  weak  shock  wave 
strikes  the  half  plane  are  determined,  in  the  linear  or 
acoustic  approximation.   The  problem  is  solved  by  means  of 
an  extension  of  the  method  of  conical  fields  which  is  well 
known  in  the  linearized  theory  of  supersonic  aerodynamics. 
In  this  paper  we  treat  the  case  where  the  undistvirbed 
magnetic  field  is  parallel  to  the  half-plane  obstacle. 

The  results  show  the  development  in  time,  toward  the 
ultimate  steady  state,  of  the  phenomena  associated  with  the 
reflection  of  a  wave  from  a  full  plane.   Ii  this  way  a 
peculiar  feature  of  hydromagnetic  waves  is  elucidated. 
Namely,  that  sometimes  the  reflected  wave  is  on  what  appears 
to  be  the  wrong  side  of  the  obstacle,  the  opposite  side  from 
the  incident  wave. 
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NYO-93514. 
DIFFRACTION  OP  A  HYDROMAGNETIC  WAVE  BY  A  HALF  PLANE 

1.   Introduction. 

In  this  paper  we  solve  for  the  diffracted  wave  which 
results  when  a  weak  hydromagnetlc  shock  impinges  on  a 
rigid  perfectly  conducting  half  plane. 

We  assume  a  medium  which  is  a  compressible  perfect 
fluid  (no  heat  conduction,  no  viscosity)  and  a  perfect 
electrical  conductor.   In  the  undisturbed  state  the  fluid 
is  at  rest  and  is  permeated  by  a  constant  magnetic  field 
B^.   The  fluid  extends  through  all  space  around  an 
obstacle,  in  the  shape  of  a  half  plane,  made  of  a  rigid 
perfectly  conducting  material.   Then,  possibly  due  to  a 
distant  distiirbance,  a  weak  plane  discontinuous  wave  or 
weak  shock  propagates  through  the  fluid.   We  are  interested 
in  the  diffracted  wave  caused  by  the  interference  of  the 
half  plane  obstacle  after  the  incident  wave  arrives.   We 
assume  that  the  amplitude  of  the  disturbance  of  the  fluid 
is  small,  so  that  it  is  appropriate  to  employ  the  linearized 
equations  of  motion. 

We  consider  only  a  special  case  of  the  diffraction 
problem^   If  cartesian  coordinates  are  so  chosen  that  the 
obstacle  is  the  half  plane  x  >  0,  y  =  0,  then  the  case 
we  consider  can  be  defined  by  saying  that  B   is  in  the 
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direction  of  the  x-axis,  and  the  direction  of  propagation 
of  the  incident  wave  is  in  the  x  y-plane*   Other  cases  can 
be  treated  by  our  method.   This  treatment,  however,  is 
considerably  more  complicated,  particularly  if  B   is  not 
parallel  to  the  obstacle. 

Oiir  method  is  based  on  the  fact  that  considerations 
of  symmetry  and  similarity  show  that  the  solution  of  the 
problem  depends  only  on  two  cartesian  coordinates  x,  y 
and  the  time  t,  and  only  by  way  of  the  ratios  x/t,  y/t» 
The  field  of  flow  is  a  conical  field  in  (x,y,t)  space. 
We  introduce  formulas  which  represent  the  solution  in 
terms  of  a  single  analytic  function  of  a  complex  variable 
z.  Here  the  z-plane  is  obtained  from  the  plane  of 
(x/t,  y/t)  by  a  non-conformal  mapping.   Using  this  repre- 
sentation, we  obtain  a  problem  of  potential  theory  in  the 
z-plane.   We  solve  the  potential  problem  and  thereby  the 
original  wave  diffraction  problem. 

The  method  of  conical  fields  has  been  used  to  solve 
a  great  variety  of  problems  in  the  linearized  theory  of 
supersonic  flow.   This  method  was  introduced  by  Busemann   , 
We  have  adapted  this  method  to  the  linear  theory  of  hydro- 
magnetic  waves,  in  which  a  more  complicated  system  of 

differential  equations  is  involved.   Our  approach  is  more 

(2) 
closely  related  to  that  of  Goldstein  and  Ward    than  to 
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that  of  Busemanno   The  starting  point  is  the  representation 
of  an  arbitrary  solution  as  a  superposition  of  plane-wave 
solutions  (including  both  complex  and  real  plane  waves). 
The  conical-field  method  has  been  applied  to  diffraction 

problems  in  electromagnetic  theory  and  acoustics  by  Keller 

(3) 
and  Blank    «   The  present  work  is  thus  a  partial  general- 
ization of  the  worK  of  Keller  and  Blank.   In  devising  the 
correct  formulatio'^  of  the  conical-field  method  for  hydro- 
magnetic  waves,  the  author  was  guided  by  the  recent  results 
of  H,  Weltzner  ^  v   Weitzner  has  determined  explicit  formulas 
for  the  fundamental  solution  of  the  hydromagnetic  equations 
in  the  two-dimensional  case*   Examination  of  the  structure 
of  Weitzner' s  formulas  revealed  the  right  way  to  generalize 
the  conical-field  method  to  hydromagnetic-wave  problems. 
The  author  wishes  to  thank  Dr,  Ho  Weitzner  for  many 
lllxaminatlng  discussions. 
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2o   Equations  of  Motion. 

The  state  of  motion  of  the  kind  of  fluid  we  consider 
Is  specified  by 

the  pressure  p  and  density  p 

the  fluid  velocity  u 

the  magnetic  vector  B 

the  electric  vector  E 

the  electric  current  density  J 
These  quantities  are  governed  by  the  Lundqulst  eauatlons  "^  , 
as  follows: 
Equation  of  continuity: 


II  +  V  .  (p  u)  =  0  (1) 


Ponderomotive  equations: 


1^  +  (u.V)  u  +  ^Vp  +  ^  Bx  J  =  0  (2) 


Maxwell's   equations:      (neglecting  displacement    c\irrent) 


V-B  =   0  (3) 


||  +  VxE=0  (1^) 
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^■^t  =  ^J  (5) 

Ohm's  Law  for  a  perfect  conductor: 

E  +  u  X  B  =  0  (6) 

Adlabatlc  equation  of  state  for  the  fluid: 

P  =  P(p)  (7) 

We  are  concerned  with  small  disturbances  away  from  a 
state  of  rest;  so  it  is  appropriate  to  linearize  the 
Lundquist  equations  by  setting 


p  =  p^(l+p^)  J  u  =  u^   ;  B  =  B^(e+B^)         (8) 

where  p  ,  B  e  =  B   are  the  undisturbed  values  of  density 
o'   o      o  "^ 

and  magnetic  field,  and   e   is  a  unit  vector.   We  insert  (8) 
into  the  Lundquist  equations  and  neglect  terms  of  second 
degree  in  p,,  u,,  B,  and  their  derivatives.   Using  (5),  (6), 
(7)  to  eliminate  J,  E,  p,  we  obtain  the  following:   (to  avoid 
clutter  the  subscripts  are  omitted). 
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II  +  V.  u  =  0  (9) 


II  +   a^  V  p   +  A^e  v(v^B)    =   0  (lO) 


§  +  ^x(e  xu)    =   0  (11) 


V  B  =    0  (12) 


Here      a,    A      are    the   iindistiirbed  values   of  the    soiind   speed 
and  AlfvSn   speed,    respectively,    defined  by 


2   =   d£        .2        _o Q3) 


3o   Plane  Waves. 

A  plane  wave  is  a  solution  of  (9),  (10),  (11),  (12) 
in  which  p,u,B  depend  only  on  the  single  parameter  o^=  Vt  -k-x. 
Here   V  is  the  propagation  speed,  and  k   is  the  unit  vector 
in  the  direction  of  propagation.   The  properties  of  the 
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various  types  of  plane  waves  have  been  discussed  by  several 
authors,  such  as  Grad^  ' ,    and  Bazer  and  Pleischman^''' \   We 
shall  review  here,  in  a  form  adapted  to  our  purposes,  some 
of  the  pertinent  resultSo 

To  study  the  plane-wave  solutions,  one  inserts 
P  =  p(cr),  u  =  u(5),  B  =  B{6)    into  (9),  (10),  (11),  (12), 
One  obtains  a  system  of  linear  homogeneous  equations  for 
the  derivatives  p',u',B«,  with  coefficients  depending  on 
V  and  ko.  As  a  compatibility  condition  one  obtains  a 
relation  expressing,  how  V  depends  on  ko   One  also  obtains 
relations  expressing  how  the  ratios  between  p«  and  the 
various  components  of  u'  and  B'  depend  on  k. 

To  state  the  results,  it  is  convenient  to  introduce 
a  cartesian  coordinate  system.   Let  the  x-axis  be  in  the 
direction  of  B^  =  B^e,  and  let  the  y-axis  be  in  the  plane 
defined  by  B^  and  k.   Let  9  denote  the  angle  between 
k  and  e,  so  that  the  x,y,  and  z  components  of  k  are 
(cos  6,  sin  9,  0)c   Let  us  denote  the  x,y,  and  z  components 
of  B  by  B^,B  ,B^;  and  the  x,y,  and  z  components  of  u  by 
(u,v,w)o. 

The  dependence  of  V  on  9  is  given  by  the  equation 


(V^-  A^  003^9)  [V^-'{A^+e.^)V^+  A^a^cos^Q]    =   0 


(11+) 
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2 
In  general  (if  cos  9  is  not  0  or  1),  this  equation  has 

three  positive  roots  V-,Vp,V->,  with  V,>V_>V-.,  corresponding 

to  three  types  of  plane  wave  which  can  propagate  in  the 

direction  9.   The  speeds  V^  and  V~.   are  the  roots  of  the 

second  factor  of  (ll|)  and  the  corresponding  types  of  plane 

wave  are  called  the  "fast  wave"  and  the  "slow  wave".   The 

root  of  the  first  factor  is  the  speed  of  the  "Alfv^n  wave" 

or  "intermediate  wave"* 

The  dependence  of  V  on  9  may  be  illustrated  graphically 

by  the  diagram  of  Pigo  lo   (This  diagram  is  not  drawn  to 

scale,  but  is  distorted  for  convenienceo )   In  a  given 

direction  6  the  fast  wave  travels  with  speed  V.  =  OP,  where 

P,  T^  is  normal  to  OP,  and  a  tangent  to  the  outer  branch 

ABCD  of  the  curve  of  Pigo  1,  and  OP,  makes  the  angle  9  with 

the  x-ajcis  or  B  o.   The  slow  wave  moves  with  speed  V^  =  OP-,, 

o  d  ^' 

where  P-^T^  is  tangent  to  one  of  the  three-cusped  inner 
branches  of  the  curve  of  Pigo  1,  and  normal  to  OP^j  and 
the  Alfven  wave  travels  with  speed  Vp  =  OPp,  where  APp  is 
normal  to  OPp,   The  distance  OA  is  the  Alfven  speed  A, 
the  distance  Oa  is  the  sound  speed  a,  and  the  distance  Oc  is 


c  =-^ 


iKT" 
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The  diagram  of  Pigo  1  is  drawn  for  the  case  A>ao   For  the 
case  A<a  the  pict\ire  is  similar,  with  A,  a   interchanged. 
For  convenience  we  shall  always  assume  A>ao 

Consider  the  totality  of  plane  waves  which  can  propagate 
in  all  directions  in  the  x  y-plane,  and  suppose  that  the 
front  of  each  such  wave  passes  through  0  at  time  t  =  0« 
Then  at  the  time  t  =  1  the  various  wave  fronts  form  three 
families  of  straight  lineso   The  envelope  of  the  fast-wave 
fronts  is  the  curve  ABCDo   The  envelope  of  the  slow-wave 
fronts  is  made  up  of  the  two  curves  CPaE  and  GHI,  together 
with  the   straight  segment  Gc»   The  envelope  of  the 
Alfven  wave  fronts  is  the  straight  segment  CAo 

The  structure  of  the  plane  waves  of  the  various  types 
is  given  by  the  way  that  p  and  the  components  of  u,  B  vary 
through  the  wave.   It  can  be  shown  that  for  Alfven  waves 


dp  =  du  =  dv  =  dB  =  dB   =  0  (l5) 

X     y 


and 


dw:  dB  =  -Vs  cos  9  (l6) 

z 


For  a  fast  or  slow  wave  we  have 


dw  =  dB  =  0  (17) 

z 
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and 

dp:  du:  dv:  dB  :  dB  =  ^xf^Ail.  (V^-a^)cos  9.  ^.  ^^^  q._^^^   q   (^q) 


In  the  limit  where  p,u,v,  etc,  vary  abruptly  and  dls- 
continuously  from  one  set  of  values  for  CT  <  0  to  a  different 
set  of  values  for  (5"  >  0  we  have  a  plane  discontinuous  wave 
or  weak  shock.   Then  formulas  like  (16)  or  (18)  relate  the 
jumps  in  p,u,v,  etc,  across  the  wave  front. 


L|.,   The  Diffraction  Problem. 

We  now  describe  the  problem  which  will  be  solved  in 
the  present  paper.   We  assume  the  following  situation:  a 
rigid,  perfectly  conducting  half  plane,  which  we  take  to 
be  the  plane  y  =  0,  x  >  0  is  immersed  In  a  space  otherwise 
full  of  a  perfectly  conducting  perfect  fluid.   Initially 
the  fluid  is  at  rest  and  permated  by  a  constant  uniform 
magnetic  field  in  the  direction  of  the  x-axls.   Then  a 
plane  wave,  carrying  perturbations  of  the  density,  velocity 
and  magnetic  field,  possibly  caused  by  some  distant  dis- 
turbance, propagates  toward  the  half  plane  in  a  direction 
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normal  to  the  edge  of  the  half  plane.   This  situation  is 
illustrated  by  Fig.  2o   We  suppose  that  the  incident  wave 
is  a  weak  hydromagnetic  shock  wave.   That  is,  we  suppose 
that  the  functions  p(<5')j,  u(tf'),  B(<r)  are  step  functions, 
which  change  abruptly  from  zero  for  o"  <  0  (ahead  of  the 
wave  front)  to  constant  final  values  for  o* >  0  (behind 
the  wave  front )o   A  full  discussion  of  such  weak  dis- 
continuous hydromagnetic  waves  has  been  given  by  Bazer  and 
Pleischman  (loco  cito)o 

To  solve  our  problem  we  must  find  a  solution  of  the 
linearized  Lundquist  Equations  (9)-(l2)  which  satisfies 
the  above-described  initial  conditions  and  satisfies  the 
•appropriate  boundary  conditions  at  the  half  plane.   At  a 
rigid  perfectly  conducting  boundary,  the  normal  component 
of  the  fluid  velocity  must  be  zero,  and  the  tangential 
component  of  the  electric  field  must  be  zero.   Let  n  be 
a  unit  vector  normal  to  such  a  boTindary.   Using  Ohm's  Law 
(Equation  (6))  we  see  that  for  the  linearized  problem  these 
boundary  conditions  are 


n  .  u  =  0  (19 ) 


nK(e-Ku)=(n'u)e='(n'e)u=0  (20) 


Ik  - 


If  the  iindlsturbed  magnetic  field  is  not  parallel  to  the 
boundary,  so  that 


n  .  e  /t  0 


then  (19)  and  (20)  are  equivalent  to 


u  =  0       at  the  boundary. 

But  for  the  problem  we  are  now  considering,  the  undisturbed 
magnetic  field  is  parallel  to  the  boundary,  so  that 

n  »  e  =  0 

Thus  for  our  problem  the  boundary  conditions  reduce  to 

n.u=v  =  0     aty  =  0,  x>0  (21) 

By  (1^),  (16),  we  see  that  if  the  incident  wave  is 
an  Alfven  wave  it  will  not  be  affected  at  all  by  the  half 
plane o   Thus  we  need  consider  only  incident  fast  waves 
and  slow  waves©   Then  initially  w  and  B   are  zero.   Also, 

2j 
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the  other  quantities  characterizing  the  flow  do  not  depend 
on  the  third  cartesian  coordinate  z.   So  we  assume  that 
w  and  B^  are  always  zero  and  that  p,u,v,B  ,B  depend  only 
o^  x,y,t   which  assumption  is  consistent  with  the  various 
conditions  of  the  problem.   Furthermore,  if  we  measure 
time  from  the  Instant  at  which  the  incident  wave  strikes 
the  edge  0  of  the  half  plane,  then  it  is  self-consistent 
to  assume  that  p,u,v,B^,B  depend  only  on  x/t,  y/tj 
so  that  snapshots  of  the  flow  at  successive  instants 
would  be  one  picture,  being  enlarged  uniformly  at  a 
uniform  rate. 

We  now  indicate  the  situation  for  t  >  0.   Pig.  3 
pertains  to  the  case  where  the  incident  wave  is  a  fast 
wave.   The  line  MT^P^KN  indicates  the  location  of  the 
incident  wave  front  were  it  not  for  the  interference  of 
the  half  plane.   Actually,  as  soon  as  the  incident  wave 
strikes  the  half  plane,  waves  are  sent  out  carrying  the 
news.   The  line  KL  is  the  front  of  the  wave  which  carries 
the  news  that  the  half  plane  is  there.   The  news  that  the 
half  plane  has  an  end,  at  0,  is  carried  by  all  possible 
waves  which  can  originate  at  0  at  time  t  =  0,   So  the 
front  of  the  region  where  this  information  is  disseminated 
is  the  curve  AT^BCDLA*   Hence  the  picture  shown  in  Fig.  3, 
In  the  region  behind  MT^BCDLKN  we  have  the  incident  wave, 
as  yet  unaffected  by  the  interference  of  the  half  plane. 
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In  the  region  KAL  we  have  a  flow  identical  with  that  which 
would  obtain  if  the  half  plane  had  no  edge  -  i.e.  the  flow 
that  would  obtain  if  the  half  plane  were  a  full  plane. 
Thus  in  KAL  the  flow  is  like  that  after  reflection  of  an 
incident  fast  wave  by  a  plane.   For  similar  reasons,  there 
is  no  disturbance  in  the  region  KT^Ao   Thus  KL  may  be 
identified  as  the  reflected  wave  front,  and  KT^  as  the 
shadow  of  the  half  plane.   Inside  AT^BCDLA  the  fluid  has 
received  the  intelligence  that  there  is  a  half  plane  and 
that  it  has  an  edge  at  0;  so  in  this  region  we  have  the 
diffracted  wave.   These  remarks  will  be  substantiated  by 
exact  calculations  which  will  follow. 

If  the  incident  wave  is  a  slow  wave  we  have  a 
different  picture,  shown  in  Fig,  I4.,   The  diffracted  wave 
in  this  case  outruns  the  incident  wave    MV  and  LN, 
Discontinuities  in  p,u,v,B  ,B   occur  across  the  segments 
VK,UK,T^L;  and  there  is  no  discontinuity  across  T^K,   Hence 
here  UK  should  be  called  the  reflected  wave,  and  T^K  the 
shadow.   The  shadow  is  on  the  same  side  of  the  obstacle 
as  the  distxarbance  giving  rise  to  the  incident  wave,  and 
the  reflected  wave  is  on  the  opposite  side  of  the  obstaclel 
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5.   Two-dimensional  Equations  of  Motion. 

To  prepare  for  a  precise  mathematical  formulation  of 
our  problem  and  for  the  calculations  which  will  solve  it, 
we  procede  with  a  more  detailed  study  of  the  eouatiois  of 
motion.   As  we  have  seen,  one  may  expect  that  for  our 
problem  we  can  assiime  that  the  components  w,B   of  u,B  are 
zero,  and  that  p,u,v,B^,B   depend  only  on  x,y,t.   Under 
these  conditions.  Equations  (9)-(12)  can  be  written  as 
follows: 


at   ax   ay  ^22; 


it  ^  ^  |^=  0  {23) 


^B 

Dt^^ff=0  (25) 


aB 

St    ax   ^  (26) 


aB    9B 

-J^  +  — 1=  0 
ax   97 


(27) 
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We  note  that  solutions  of  this  system  can  be  repre- 
sented In  terms  of  a  single  function  \|/,  a  sort  of  potential. 
By  (27)p  we  see  that  \|/  exists  such  that 


B  =  -lir   ^  (28) 


B^  =  ^V_.  (29) 

y    xyt 


(subscripts  appended  to  ■^   denote  partial  derivatives)* 

Then  from  (25),  (26)^  (21;),  (23)  we  see  that  these  equations 

are  satisfied  if 


^ytt 


(30) 


p-^hft  [^tt-^'^^xx^  V^^  ^^^^ 


-=  -Ie  t*tt-^'(*xx^  V^^  ^^^^ 


Then  from  (22)  we  obtain  the  following  equation  for  \|/: 


dS.   .(a2.  a2)^  (^  +  ^).  a^A^  ^  (^fi  ^  2%)  =  o   (33) 
at^         3t^  ax'^   ay"^       gx"^  ax^   ay"^ 


=  19  - 


We  must  make  a  few  remarks  here  about  the  propagetion 
of  discontinuities  in  a  solution  of  the  equations.   This 

subject  has  been  dealt  with  in  detail  by  Bazer  and 

(7) 
Fleischman   ,  Let  us  suppose  we  have  a  solution  of 

(22)- (27)  which  tends  to  different  limits  as  the  point 

(x,7,t)  tends  to  the  two  different  sides  of  a  surface  S  in 

(x,y,t)  space,  so  that  p,u,v,B^,B   have  jump  discontinuities 

across  S,   Of  course,    then  the  differeitial  equations  break 

dovm.   One  must  use  the  laws  of  conservation  of  mass, 

momentum,  and  energy,  and  Ampere's  law,  which  are  integral 

relations.   In  this  way  it  can  be  shown  that  at  a  discontinuity 

surface  the  jump  in  p,u,v,B^,B  ,  which  we  denote  by  [p],  [u] , 

[v],  [B^],  [B  ]  must  satisfy  the  following  conditions: 


n^[p]    +  n^[u]    +  n   [v]    =0  (31^) 


n^[u]    +   a^n^[p]    =   0  (35) 


n    [v]    +   a^n   [p]    +  A^(n^[B^]    -   n^[B^] )    =  0  (36) 


n^-tB^]    -   ny[v]    =  0  (37) 


n^[By]    +  n^[v]    =  0  (38) 


^x^^x^    ■"  ^y^Sy^    =0  (39) 
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there  n.  an  ,n  are  the  components  of  a  vector  normal  to 
t'  X*  7  ^ 

S  In  (x,y,t)  spacco   One  elimination,  it  follows  that 


n^  -    (A^  +  a2)n^(n2  +  n^)  +  A^aVCn^  +  n^)  =  0    (1|.0) 


That  is,  S  is  a  characteristic  surface  of  the  system 
(22)-(27)  and  of  the  single  equation  (33)« 

A  characteristic  surface  can  be  shown  to  be  the  envelope 
of  characteristic  planes,  which  are  fast  or  slow  plane-wave 
frontSo   A  characteristic  surface  of  special  interest  is 
the  characteristic  cone,  which  is  the  envelope  of  all 
characteristic  planes  through  a  point  in  (x,y,t)  spaceo 
This  the  characteristic  cone  with  vertex  (0,0,0)  is  repre- 
sented in  the  plane  of  (x/t,y/t)  by  the  curves  ABCD,  EOF, 
GHI  in  Pigo  I4-0 


6«   Mathematical  Formulationo,   Uniqueness  Theorem. 

We  formulate  the  boundary  problem  for  the  reflected 
and  diffracted  wavesj  that  is,  for  the  total  solution 
minus  the  incident  wavco   First  we  state  regularity  con- 
ditions required  of  p,u,v,B  ,B  ,   We  assume  that  these 
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functions  are  differentiable  in  the  open  region  R  in 
(x,y,t)  space  gotten  by  subtracting  certain  characteristic 
surfaces  and  the  half  plane  x  >  0,  y  =  0  from  all  of 
(xyt)  space.   We  assume  that  the  functions  p,u,  etc,  tend 
to  continuous  limits  as  (x,y,t)  approaches  either  side  of 
one  of  these  surfaces,  except  at  certain  lines  L . ,   In 
Fig.  3  the  discontinuity  surfaces  are:  the  characteristic 
coneABCD,  CFE,  GHI,  the  plane  wave  fronts  KT,,  KL,  the 
half  plane  OK  and  in  addition  the  plane  segment  OG  —  which 
is  characteristic.   The  lines  L.  are  represented  by 
K,T, ,L,G,  and  0,  which  last  is  the  edge  of  the  half  plane. 
We  assume  that  near  the  lines  L.,p,u,v,  etc,  are  bounded 
—  except  near  the  line  0,   Near  0  we  assume  p,u,  and  v  are 
bounded,  but  that  B   and  B  may  become  infinite  in  such  a 

way  that  for  any  t^  >  0,  |B  |<  C{t    )/i~r,     |B  |<  G(t  )//?" 

/~2      2 
for  t  <  t  ,  where  r  = /x  +y  , 
o 

We  assume  that  in  R  the  equations  (22)-(27)  are 
satisfied.   We  assume  that  at  a  characteristic  surface  of 
discontinuity  the  relations  (3l;)-(39)  are  satisfied. 

We  assume  that  v  satisfies  on  x  >  0  y  =  0  the  following 
boundary  condition: 

for  y=iO,      0<x<  V/cos   0,    v  =   1 

for  y  =   i  0,      V/cos   0   <  X,    V  =   0 
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Here  of  course,  we  assvime  that  the  incident  wave  travels 
with  speed  V  in  direction  9,  and  carries  a  value 

V  =  -  1. 


We  also  assume,  throughout  that  0  <  9  <  90  .   This  is 


purely  a  matter  of  convenience. 
Finally,  we  assume  that 


for  t<0,  p  =  u=v  =  w  =  B=B=0, 

■X.    y 


We  also  assume  that  for  t  =  t^  >  0  we  have  p  =  u=v=w=B^=B=  0 
outside  of  some  bo\inded  region  of  the  xy  plane. 

This  last  assumption,  which  says  merely  that  hydromagnetic 
waves  travel  with  finite  speed,  can  be  proved  to  be  a  con- 
seauence  of  the  other  assumptions,  but  we  shall  not  stop 
here  to  give  the  proof. 

We  now  prove  the  uniqueness  theorem;  that  the  above 
problem  can  not  have  more  than  one  solution.   It  is  sufficient 
to  show  that  the  difference  of  two  solutions  is  zero;  and 
of  course  the  difference  of  two  solutions  satisfies  the 
homogeneous  boundary  condition 

for  y  =  ±  0,  X  >  0,  V  =  0  (l|l) 
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The  proof  follows  from  conservation  of  energy.   We 
multiply  (22)  by  a^p,  (23)  and  (21+)  by  u  and  v,  (25)  and 
(26)  by  A  B^  and  A  B  ,  and  add  the  results.   We  obtain 
a  relation  which  can  be  written  as  follows: 


^[1,V^1u2./,.1,2,b2.b2,, 


+  ^  [a^pu  -  A^vB  1  +  ^[v(a^p  +  A^B^)]  =  0     (1^2) 


We  now  integrate  this  over  all  of  (xyt)  space  with  t  <  t  . 

We  obtain  boundary  terms  at  each  surface  S^  where  p.u.v.B  .B 

i        '  *  *  X*  y 

do  not  have  continuous  derivatives,  or  are  discontinuous* 
These  boundary  terms  arise  from  the  second  two  terms  of 
(I4.2),   The  singularities  at  L.  contribute  nothing,  by  virtue 
of  the  assiomptions  which  were  made  about  the  behavior  of 
p,u,v,B  ,B   near  L.,   The  boundary  term  arising  from  the 
half  plane  x  >  0,  y  =  0  is  zero  by  virtue  of  the  homogeneous 
boundary  condition  (1+1), 

Also,  the  boundary  term  from  a  characteristic  surface 
S  is  zero.   For  the  integrand  of  the  surface  integral 
arising  at  S  is 
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njl  a^p^   +  |(u2  +  7^)    +1  A^d^   +  B^)] 


+  n    [a^pu  -  A^vB   ]    +  n    [v(a^p  +  A^B    ) ] 


-  n^[|  a^p^  +  I   (u^   +   v^)    +  I  a2(B^  +  B^) ] 


-  n    [a^pu  -   A^vB    ]    -  n    [v(a^p   +  A^B^)]  (W) 

X      "^  y         y  ^ 


where  (n  ,n  ,n. )  are  the  components  of  the  unit  normal  to 
X*  y'  t 

S  in  (x,y,t)  space  and  p,  p;  u,  u  etc.  are  the  values  to 
which  p,  u  etc.  tend  on  the  two  sides  of  S,   And,  the 
expression  (14-3)  can  be  rewritten  as  follows: 


I  a2(p   -H   p; 


)    n^[p]    +  n^[u]    +  ny[v]j 


+  i   (u  +  u) 


r  2 

n.  [u]    H-  a     n^tP^ 


+  I   (v  +   vjfn^Lv]    +   a^n^Lp]    +  A^n^lB^]    -  A^m^[By]j 


+  i   (B^-^B^)I 


n,[B^]    ^  ny[v]] 


+  ^   (B     +  B    )    n^[B    ] 


2       y         y  I     t^   7 


[n^[B^]    -   n^[v]] 
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where  [p],  [u]  are  the  Jumps  "p  -  p,  u  -  u  etc.  in  p,  u,  etc. 
Thus  we  see,  by  (3l;)-(39)  that  this  boiindary  terra  is  zero. 

Hence  after  the  integration  has  been  carried  out,  we 
are  left  with 


]\ 


[|  a^p^  +  |(u^  +  v^)+  I  A^(B^  +  B^)]dxdy  =  0      (Ul|) 


and  (I4I1.)  implies 


p  =  u=v=B  =B  =0 
X    y 


which  proves  the  uniqueness  theorem, 


7»   Representation  of  a  Conical  Field  by  Plane  Waves, 

We  now  develop  the  mathematical  machinery  by  which  our 
problem  will  be  solved. 

Let  us  consider  a  conical  field;  that  is,  a  solution 
of  (22)-(27)  in  which  p,u,v,B^,B   depend  only  on  x/t,  y/t, 
or,  put  another  way,  in  which  p,u,v,B  ,B   are  homogeneous 
functions  of  x,y,t,  of  degree  zero.   We  can  construct  such 
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a  conical  field  by  superposing  plane  waves  which  are  conical. 
Indeed,  such  a  solution  is  eriven  by        ■      * 


(t,x,y)  =  -  ^  J  tt  -  :^(3)x  -  ti(s)y]^  f(s)ds  (1^5) 


vhere  A(s),    ^.(s)    are   functions   of   s   related   by 


F{-K,[l)   =   1   -    (a^  +  A^){\^  +   [i^)+   R^kh^U^  +  [i^)=   0        (l|6) 


and  where  z  is  a  function  of  (x/t),  (y/t)  such  that 

t  -  A(z)x  -  n(z)y  =  0  .        .         (1^7) 

The  function  f   defined  by  iliS) ,    ik^) ,    ihl)    satisfies  Equst  ion 
(33)s  as  can  be  verified  formally  by  direct  differentiation, 
making  use  of  the  following  result  of  differentiating  ikl)' 


St       dx       ay 


The  solution  p,u,v,B  ,B   is  given  by  third  derivatives  of 
\|/  according  to  Equations  (28)-(32)o   Since  \|/  evidently  is 
homogeneous  in  x,y,t  of  degree  3,  the  third  derivatives 
p,u,v,B  ,B   are  homogeneous  of  degree  zero, 

^  y 
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Equations  (1^.6)  and  {1^7)    say  that  l/X,  l/ti  are  the 
intercepts  with  the  x/tp  y/t  axes  of  a  tangent  through 
P  E  (x/t,  y/t)  to  the  characteristic  cone  of  Pig,  1,   In 
general  through  a  given  point  P  =  (x/t,  y/t)  there  are 
four  such  tangentso   (Eliminating  \i   from  (I4.6),  (l;7)  one 
obtains  a  quartic  equation  for  \,)      One  can  see  from  Fig,  1 
that  for  P  outside  the  outer  branch  of  the  characteristic 
cone  the  four  tangents  are  real,  being  two  fast  waves  and 
two  slow  waves o   For  P  inside  the  inner  branch,  the  four 
tangents  are  real,  being  four  slow  waveso   For  P  between 
the  inner  and  outer  branch,  one  has  only  two  real  tangents  - 
slow  waveSo   The  other  tangents  are  complexo   (This  is,  \ 
and  ^L  are  complexo ) 

Corresponding  to  the  foior  different  tangents,  there 
are  four  possible  formulas  (l+S)©   We  may  surmise  that  by 
adding  formulas  of  the  foiir  types  one  obtains  the  general 
conical  field,  expressed  in  terms  of  four  functions  f^ 
This  is  true,  although  we  shall  not  prove  it  here.   To 
solve  our  problem  it  suffices  to  consider  only  the  two 
complex  waves  associated  with  a  point  between  the  outer 
and  inner  branches  of  the  characteristic  cone.   We  may 
remark  in  passing  that  for  a  conical-field  problem  in 
which  part  of  the  boundary  is  not  parallel  to  the  undis- 
tijrbed  magnetic  field,  the  waves  corresponding  to  the  real 
tangents  also  occur,  and  must  be  consideredo 
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We  now  writ©  dowi  a  formula  for  conical  fields  which 
represent  a  solution  at  points  through  which  two  tangents 
are  compleXo   We  define  X(z)  and  p,(z)  in  the  region  in  the 
z-plane  outside  the  slits  ca,  AB,  GK,  CB  (see  Pig.  5) 

by 

■k  =   l/z 

(i|8) 


^=^7K^^-.^/ 


(z  -A  ) (z  -a  ) 

(    2   2. 
(z  -c  ) 


The  square  root  occvu?lng  in  \i   is  to  be  defined  so  as  to 
be  positive  real  for  z  real  and  on  the  upper  part  of  the 
slit  from  z  =  A  to  z  =  +oo  • 

In  (l;8),  c  is  defined  by 

o  2,2  • 

«2  _  a  A 
c  -  -^ — ^  . 

a  +A 

The  fionctions  {\\.Q)    satisfy  the  relation  (i|6),  as  can 
be  verified  directly. 

The  relation  between  (x/t,  y/t)  and  z  is  then  given 
by  ikl),    or 

(i)(Vt)*(|j)(i)/2lA!ru!IZ)  (y/t,  =  1.    (1^9) 
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The  real  and  imaginary  parts  of  Equation  (14-9)  are 
two  equations  for  (x/t),  (y/t)  in  terms  of  z.   This 
equation  defines  a  dif ferentiable  one-to-one  mapping  of 
the  open  region  R-.  between  the  outer  branch  ABCD  and 
inner  branches  cPaE,  GHKI  of  the  characteristic  cone  of 
Pigo  1  into  the  slit  z-plane  of  Fig,  5«   (We  shall  prove 
this  in  the  AppendlXo)   In  this  mapping  the  arc  BAD  in 
Fig.  1  is  mapped  into  the  slit  BAD  in  Fig.  5  (where  B,D 
are  mapped  into  the  point  at  infinity  in  the  z-plane.) 
The  arc  BCD  is  mapped  into  the  slit  BCD  in  the  z-plane. 
The  c\arve  cEaF  is  mapped  into  the  Flit  cEaF,  and  the 
curve  GHKI  is  mapped  into  the  slit  GHKI,   A  point  M  on 
the  axis  in  R,  (see  Fig,  1)  is  mapped  into  a  point  M 
on  the  real  line  in  Fig,  6,   The  origin  0  in  Fig,  1  is 
mapped  into  the  origin  0  in  Fig.  5o   Also  the  distance 
OM  is  not  changed  by  the  mapping.   That  is,  the  distance 
OM  in  Figure  1  is  equal  to  the  distance  OM  in  Fig,  ^, 
Also,  near  0  the  mapping  of  the  (x/t),  (y/t)  plane  Into 
the  z-plane  approximates  the  identity  mapping,  in  the  sense 
that  (if  z  =  ^  +  i"^)  we  have 

C  =  (x/t)[l+0(lzl'+)] 

(50) 
7)=   (y/t)[l+o(U|^)] 
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The  Intuitive  reason  for  all  this  is  that  z  has  the 
following  geometric  interpretation?  to  determine  z,  draw 
through  P  E.  (x/ts,y/t)  one  of  the  complex  tangents  to  the 
characteristic  cone.   Then  z  is  the  intercept  of  that 
tangent  on  the  x/t  -  axiSo   Since  the  point  M  is  already 
on  the  x/t  axis,  any  line  through  M  intercepts  the  x/t 
axis  at  M.   Hence  OM  in  Fig„  5  is  equal  to  OM  in  Pigo  lo 
As  the  point  P  approaches  the  outer  branch  ABCD  of  the 
characteristic  cone  at  T-,  ,  the  two  imaginary  tangents 
become  equal  and  therefore  realo   Hence  T^  is  mapped  on 
to  the  foot  U^  of  the  tangent  T^U^o   Similarly  if  P 
approaches  T^  on  an  inner  branch  aEcF  of  the  characteristic 
cone,  its  image  in  the  z-plane  tends  to  a  point  TJ-,  on  the 
real  axis  with  OU^  (Figo  5)  =  OU^  (Figo  1)  where  T^U^  is 
tangent  to  aEcFo   All  these  statments  about  the  properties 
of  the  mapping  will  be  proved  in  the  AppendiXo 

Now  we  are  prepared  to  state  the  formula  with  which 
we  will  solve  the  boundary  problem.   This  formula  is 


oo 
\|>(x,y,t)  =  -|  ImJ  [t  -  \(s)x  -  M.(s)y]-^f(s)ds    (51) 

z 

where 


(  z  <=c  ) 
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and  z  is  related  to  x/t,  y/t  by 

X(z)(x/t)  +  ti(z)(y/t)  =  1  (53) 

Here  f(z)  is  required  to  be  single-valued  in  the  z-plane 
slit  as  follows:   there  are  slits  BC,  KG  on  the  negative 
real  axis,  from  -oo  to  -A  and  from  -a  to  -Co   Also  there 
is  a  slit  from  0  to  oo  on  the  positive  real  axiSe   (The 
positive  real  axis  is  in  effect,  the  image  of  the  half 
plane  obstacle). 

The  path  of  integration  in  (50)  is  arbitrary  except 
that  we  assume  the  path  can  cross  the  real  ajcis  only  on 
the  segment  CK,  between  z  =  -A  and  z  =  -a©   To  ensiore 
that  ^    so  defined  is  Independent  of  the  path,  we  require 
that 


f(2)  =  0(l/z^)  as  z  -*  00  iSk) 


f(z)  is  real  and  integrable  on  the  slit  BCD     (55) 


To  complete  the  solution,  we  must  now  choose  f(z) 
so  that  p,u,v,B  ,B   as  obtained  from  (50)  by  (28)-(32) 
satisfy  the  regularity  conditions,  and  so  that  v  satisfies 
the  boundary  conditionso 
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8.      Solution  of   the    Problenio 

In  this  section  we  recapitulate  the  formulation  of 
the  diffraction  problem,  and  then  propose  and  verify 
the  solution.   We  state  the  conditions  of  the  problem  in 
a  stronger  form  than  heretofore »   We  shall  assume  that 
the  angle  of  incidence  0  (see  Figo  2)  is  between  0  and 
90°  and  that  A  >  a,  though  it  is  a  simple  matter  to  give 
the  solution  in  other  cases. 

Suppose  the  incident  wave  is  a  fast  wave.   Then  the 

problem  is  to  find  quantities  p,u,v,B  ,B   depending  only 

X  y 

on  (x/t,  y/t)  which  satisfy  the  following  conditions: 

(1)   Except  at  the  curves  ABCD,  cEaF,  GHI,  and  the 
line  segments  OG,  OK,  KT,,  KL  (see  Figure  3)  the  quanti- 
ties are  dif ferentiable  and  satisfy  Equations  (22)-(27). 
Here  K  is  the  point 


x/t  =  V/cos  9  =  z 

o 


on  the  half  plane  obstacle  and  the  front  of  the  incident 
wave, 

(2)   Outside  the  region  bounded  by  the  two  char- 
acteristic lines  KT^p  KL  and  the  arc  T-BCDL  of  the  char- 
acteristic cone,  the  quantities  are  all  zero. 
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(3)   On  the  characteristic  cone  the  quantities  are 
continuous,  except  at  T^,L,G,  at  which  point  the  quantities 
are  bounded 

(I4.)   Across  the  line  OG,  KT,  KL  the  discontinuities 
in  the  various  quantities  satisfy  the  conditions  (31;)-(39), 
These  conditions  are  equivalent  to  the  following: 

(a)  Across  OG,  a^[p]  +  A^[B  ]  =  0,  [v]=  0,  [B  )=  0 

(b)  Across  KT,, 

[p]:  [u]:  [v]:  [B^] :  [B  ]  = 


where 


1        c    1  /(2^-A^)(z^a2) 
o  0/    (Zq-c  ) 


(Using  the  positive  square  root) 

(c)   Across  KL,  a  like  condition  holds,  with 


,.1         ,c  ,  lxPo-A'>'zg-a2) 

'^  -  z '  1^  -  -  '3>  I V — rz"^^ — 

o  01    (z  -c  ) 
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(5)   At  the  half  plane,  on  OK,  the  ouantlties  tend 
to  lirnits  for  y/t  -^  +  0  or  y/t  — ^  -  0,  and 

for  0  <  x/t  <  z 

o 

we  have  - 

V  =---  1 

At  one  point  0  there  is  a  singularity.   Near  0,  p,u,v, 
are  bounded,  but 


Sx  =  0| 


it 


and  B  =  0 


^?^'       ' 


/Vt ^\ 


V 


Ix   +y 


(6)   In  each  of  the  regions  AKT^,  AKL,  aEc,  aPc,  GHI, 

the  quantities  p,u,v,B  ,B   are  constantSc   (However,  the 

X  y 

constants  are  different  for  the  different  regions.) 

With  these  assumptions  one  has  boundary  conditions 
on  the  inner  and  outer  boundaries  of  R,  ,  and  corresponding 
boundary  conditions  in  the  z-plane,  which  we  now  state, 
(See  Pig.  5) 


(7)   on  both  sides  of  BCD, 


p  =  u=v=B  =B  =0 
X     y 
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(8)  on  both  sides  of  KG, 

p  =  constant,  u  =  constant,  v  =  constant, 
B  =  constant,  B  =  constant 

-^  y 

(9)  Across   OG,    v   and  B      are    continuous    and  the    jumps 

V 

in  p  and  B   are  related  by 

a^[p]  +  A^[B  ]  =  0 


(10)   On^=  ±  0,  0  <  C  <  z^  ,  we  have 

V  =  1 
On 77=  *  0,  ^  >  z  we  have 


p  =  0,  u=0,  v  =  0,  B   =0,  B   =0 


(11)   Onr^=  +  0,  c  <  C  <  a,  the  quantities  p,u,B  ,B  ,v 
are  constant. 

On  79=  -  0,  c  <  ■?<a,  the  quantities  p,u,B  ,B  ,v  are 
constant  (but  not  necessarily  equal  to  their  values  for 

\=    +  0)o 

(12)  The  quantities  p,u,v,B  ,B    are  continuous 

^  y 

except  at  z  =  -c,  z  =  0,  z  =  z  ,   At  z  =  -c  and  z  =  z 
these  quantities  are  bounded.   At  z  =  0,  the  quantities 
p,u,v  are  bounded,  but 


B   =  0  '    1   \   o   _  .  /   1 
X 


TT^] '  \ -' °  [rim 
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Prom  (5l)  and  the  equations  (28)  to  (32)  expressing 
p,u,v,B  ,B   in  terms  of  third  derivatives  of  \|fo  and  using 
(52)  and  (1^.6)  we  find  that  a  possible  solution  for 
p,u,v,B  ,B   is  represented  by 


p  =  Im 


P(s)ds    ;   ||  =  Im  P(z) 


00 


f. 


u  =  Im   U(s)ds 


J   ||  =  Im  U(z) 


00 


where 


1 


V  =  Im   U(s)ds 


J   H  =  Im  V(z) 


00 


(56) 


)z  3b 

B^(s)d3  ;      -^  =  Im  B^(z) 

CD 


f^        aB 

B  =  Im   B^(s)ds   ;   -^  =  Im  B 

00 


(z) 


f<-)  =  <iA' 


C    N      Z 


>^^ 


>r~2   2w  2   2 

/(  z  •-a  ]  (  z  -c 


Y(z) 


-  37  - 


U(z)  =  {^)      ^^  -^  V(z) 


Vu^ 


a  )  ( z  -c  ) 


B,(z)  =  -  i  V(z)  . 

il  z 


We  must  now  select  a  function  V(z)  such  that  ^,  etc» 
have  the  right  behavior  on  the  real  axl3)j?=  0,  such  that 


P(z),  etc,  =  0(1/7  )  as  z  — **  00 


and  such  that  v  J\imps  from  0  to  1  as  C  crosses  the  point 


?  =  z  =  z^  . 


The  last  requirement  means  that  V{z)  has  a  simple  pole  at 

z  =  Z^o 

o 

It  can  easily  be  verified  that  all  requirements  are 


met  by  the  choice 
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v(z)=Q-y72 TTT T72 (^Q^ 


where 


z 


1/2,      ^    a/2,      ^,n1/2 
'     ( z_+c)    '     ( z.+A)    ' 


Then  we   have 


1/2,    , a/2 

P(z)    =    (f^)    Q  ,,l        ^^-J)l    (60) 

^^  (z-a)^/^(z-c)^/^(z+c)(z-z    ) 

9 


1/2 

U(z)    =    (f^)    Q  ^72 1/2'"^^    1/2 (^1) 

^  z-^/^(z-a)-^/^(z-c)^/^(z+c)(z-z^) 


B,(z)    =    (f^)    Q    (^J|^<^-^j!f(^-^)^^^  (62) 

^  ^  z^/^(z-c)^/2(^^^)(^_^^) 


/    .    a/2 
B  (z)  =-Q-T72 ttJ^"^ T72 .  (63) 

2  z^/^(z+c)^/^(z+A)-^/^(z-z    ) 

o 


Condition    (9)    on   the    Jumps    in  p,B      across   OG   Is   easily 


verified.      Also    it    is   easily  verified    that   the    Jumps   in 
P,u,v,B^,B^  J 
(I^b)    and    (3) 


p,u,v,B    ,B      at   z   =   z      are    in   accordance   with  conditions 
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Thus  p,u,v,B  B   are  obtained  (by  integrating  these 
formulas,  in  accordance  with  (56),  as  functions  of  z.   The 
functions  involved  are  elliptic  integrals  of  the  third 
kind.   To  obtain  p,u,v,B  ,B  as  functions  of  (x/t),  (y/t) 
we  must  obtain  z  from  (52),  {S3).      In  general  this  involves 
solving  a  biquadratic  equation.   However,  for  y/t  =  0, 
which  is  an  interesting  special  case,  matters  are  simpler, 
for  here,  as  has  been  stated, 

z  =  x/t  , 

If  the  incident  wave  is  slow,  the  solution  is  given 
by  the  same  formulas  (58)-(63)  except  that  for  a  slow 
incident  wave  we  have 


c  <  z   <  a  , 
o 


9.   Discussion  of  the  Results. 

We  now  mention  some  features  of  the  pattern  of  the 
reflected  and  diffracted  waves  that  have  been  revealed  by 
our  analysis. 
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If  the  incident  wave  is  fast,  we  have  the  situation 
illustrated  in  Pigo  3o 

Ahead  of  MT^AK  the  fluid  is  at  rest. 

Behind  MT^CDLKN  we  have  the  incident  wave  as  yet 
undisturbed  by  the  half  plane. 

In  KAL  the  flow  is  a  constant  state  (that  is,  p,u,v,B  ,B 

■^  y 

are  constant  throughout  KAL),   The  flow  in  KAL  is  Just  like 
the  flow  for  a  fast  wave  reflected  from  a  full  plane. 

In  GHI,  and  also  in  CaE,  and  in  caP  the  flow  is  a 
constant  stateo 

At  OG  there  is  a  contact  discontinuity.   Across  OG 
the  parallel  component  of  velocity  u  is  discontinuous,  and 
also  the  density  p  and  the  component  B  of  the  magnetic 
fieldo   Since  p  is  discontinuous,  so  is  the  pressure.   Since 
B   is  discontinuous „  a  sheet  cvirrent  flows  in  the  z  direction 
in  OGo   The  pressure  Jump  is  balanced  by  the  force  exerted 
by  B   on  this  sheet  current,  and  so  we  have  equilibriijm. 
This  is  the  content  of  condition  (l|a). 

At  the  edge  0  of  the  half  plane,  p,u,v,  are  continuous 
(on  the  upoer  side  and  on  the  lower  side),  but  B  ,B   become 
infinite  like  the  reciprocal  square  root  of  the  distance 
from  the  edgCo 

If  the  incident  wave  is  slow,  we  have  the  situation 
illustrated  by  Flgo  ho 

Ahead  of  MVBALN  the  fluid  is  at  rest. 


m   - 


Behind  MVCDLN  the  Incident  wave  is  as  yet  undlstiirbed 
by  the  half  plane. 

The  quantities  p,u,v,B^,B   suffer  jumps  across 
VK,  T^L  —  which  lines  are  the  front  of  the  incident  wave 
propagating  through  the  diffracted  wave  which  in  this 
case  outruns  the  incident  wave. 

Across  KU  the  flow  has  a  jump  discontinuity.   The 
line  KU  is  the  front  of  the  reflected  wave,  propagating 
through  the  diffracted  wave. 

The  flow  is  continuous  across  KT,,  which  is  the  shadow. 

In  the  following  regions  the  flow  Is  a  constant  state: 
GHI,  cKW,  WKU,  EUKa,  cKaPT^,  T^PF. 

A  remarkable  feature  of  the  results  is  the  occurrence 
of  regions  of  constant  state  inside  the  cusped  curves  GHI, 
c7:P.   This  Is  another  manifestation  of  a  phenomenon  dis- 
covered by  Weltzner  ^'    for  a  different  problem.   Weltzner 
found  the  fundamental  solution  of  the  system  (22)-(27), 
(That  is,  the  solution  of  (22)-(27)  modified  by  replacing 
the  zeros  on  the  right-hand  sides  by  Dlrac  delta  functions.) 
Weltzner  found  that  the  interiors  of  the  cusped  curves  are 
lacunae;  that  is,  the  fundamental  solution  is  Identically 
zero  in  these  regions. 
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Appendix.   Mapping  of  the  (x/t).  (y/t) -Plane  into  the  z-Plane. 


In  this  Appendix  we  wish  to  show  that  the  equations 


\(z)(f)  +  ^(z)(f)  =  1  (6H.) 


ll(z^-a^)(z^-A^) 


aA   zx/   (22_^2j 


define  a  dlfferentlable  one-to-one  mapping  of  the  z-plane 
silt  as  in  Pig.  5  Into  a  region  R^  of  the  (x/t),  (y/t)  plane 
(shown  in  Pig.  1),  with  the  following  properties: 

(1)  A  point  T  on  the  boundary  of  R,  is  mapped  into 
the  point  U  where  the  tangent  TU  cuts  the  (x/t)  axis, 

(2)  Points  in  R,  for  which  y/t  =  0  are  invariant. 

(3)  Near  0  the  mapping  approximates  the  identity 
mapping  in  that 

z  =  (f  +  1  f)[l  +  0(|zl^)]  .  (66) 


(I4.)      The   boiindlng   cxirves   of   R,    determine    characteristic 
siirf  aces. 
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To  show  that  the  mapping  is  one-to-one  we  show  first 
that  for  a  complex  value  of  z,  the  point  (x/t,  y/t )  is 
determined  uniquely,  and  is  a  dif f erentiable  function  of 
C, )?  where  C,>7  are  the  real  and  imaginary  parts  of  z«   Let 
X,  ,  \i^    be  the  real  parts  of  X,  \i.   and  \p,    [ip   the  imaginary 
parts.   The  real  and  imaginary  parts  of  (61|)  yield  two 
linear  equations,  and  the  solution  is 


t   X^|j.2->^2^1 


t   ^i^^2"'^2'^l 


We  must  verify  that  the  denominator  is  not  zero.   The 
denominator  is 


(66) 


2    2      ^l"^  ^^2 


which  is 


/aAw  2  ,   2-1  T  \^/    (z  -c  ) 
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The  imaginary  part  of  a  function  analytic  in  some  domain 
takes  on  its  maximum  at  the  boundary  of  the  domain.   Take 
the  domain  to  be  the  slit  z-plane  of  Pig.  5  with  small 
regions  excised  near  z  =  '^  a,  z  =  -  Ao   Near  these  points 
it  is  easily  shown  that 


^'"  iTTrfekpT  *  ^  ° 


and  on  the  slits  we  have 


r 


Im 


2   2 
z  -c 


(z  -a  ) (z  -A  ) 


^ 


=  0  . 


It  follows  that  in  the  interior  of  the  entire  slit  z-plane 
we  have 


\ 


Imv 


^  2 
z  -c 


,2      2s  I    2    .2^\ 
(z  -a  )  (z  -A  )j 


<  0  ;  so  ^^^2  "  ^2^^1  ^  °*   ^^'''^ 


Hence  a  dif ferentiable  mapping  of  the  z-plane  into  the 
(x/t),  (y/t)  plane  is  defined  by  (6I4.),  (65). 
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By  (66),  (67),  we  see  that  if  z  is  in  the  upper  half 
plane,  so  that  \^   is  negative,  then 

y/t  >  0 
and  for  z  in  the  lower  half  plane, 

y/t  <  0 

Thus  Im(z)  and  (y/t)  have  the  same  sign. 

Now  to  show  that  the  mapping  is  one-to-one,  we  show 
that  for  given  (x/t),  (y/t)  there  is  not  more  than  one 
value  of  z  in  the  interior  of  the  slit  z-plane  such  that 
(61^.),  (65)  hold. 

The  possible  such  values  of  z  are  among  the  four 
roots  of  the  biquadratic  equation 


WU){\)]^  =  [i-x(z)(|)]2  . 


This  equation  always  has  a  real  root  z  with  -a  <  z  £  -  c, 
and  a  real  root  such  that  c  <  z  <  a.   We  see  this  by  noting 
that  the  family  of  lines 


aA    z  Y 


1(2L)    ±  ^  IALZ^ZilZI^  (i)    =   1      -a  <    z   <   -c 
z^t^         aA    z  V      (^27^2T  ^t^         -L,      a  <    z   <   -c 
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covers  the  whole  (x/t),  (y/t)  plane,  since  the  intercept 
of  such  a  line  with  the  (y/t)  axis  can  have  all  values 
from  -  00  to  +  00  o   A  similar  observation  shows  the 
existence  of  a  root  with  c  <  z  <  a,  hence  the  biquadratic 
equation  can  have  only  two  roots  with 

Im  (z)  +  0, 

These  roots  are  conjxigate  complex  numbers,  and  for  only 
one  of  them  can  Im(z)  have  the  same  sign  as  (y/t).   This 
proves  the  mapping  is  one-to-one,  except  for  the  cases 
-A  <  z  <  -a,  -c<z<c,  a<z<A,   But  for  such  a  value 
of  z,  the  quantity  \i.   is  pure  imaginary,  and  X  is  real; 
that  is, 

^  =   0,    ^x^  =   0 
and  from  (66)  we  get 


7^=  0,  ^  =  z  =  (x/t) 


which  of  course  is  an  unique  value  of  z. 

By  the  way,  we  have  moved  property  (2)  above. 

Now  let  us  consider  the  mapping  of  the  boundary.   We 

consider 

-  hi   - 


with 


C  <  -A,  or   -&  <   E,  <   -c,  or  c  <  p  <  a,  or  A  <  C, 

and  assume  that  7?  la  small.   For 77=  0  the  point  (K,7/)    Is 
on  some  part  of  the  boundary  of  the  region  in  the  z-plane. 
Also,  f  or  )7  =  0,  A  and  \i   are  real.   Expanding  in  powers  of 
^  we  see  that  for  small  >?  the  real  and  imaginary  parts  of 
{6l\.)    are 


X(^)(|)+  ^(^)(y/t)  =  1 


(67) 


if^^^^^^f)  ^^^(^^(?)]  =  0  • 


Here  the  sign  of  \i   will  depend  on  whether  )?  is  positive 
or  negative  —  that  is  on  from  which  side  the  point  (?,^) 
tends  to  the  slit  as  |)7|  ^-»-0,   Prom  (6?)  we  obtain  the 
following  equations,  which  determine  the  image,  in  the 
(x/t),  (y/t)  plane,  of  the  bounding  slits  in  the  z-plane: 

X(C)(x/t)+  ti(^)(y/t)  =  1 

(68) 
\'(?)(x/t)+  ti'(?)(y/t)  =  0 
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The  determinant  of  these  equations  is  not  zeroj  for  this 
determinant  is 


3 


Xti'  -  txX'  =  7.2  ^  (^)  =  ^^   [(^)2] 


and  direct  computation  shows  that 


^  [(ixA)2]  =  0 


only  for 


2 

a' 


la  + 


and,  since  A  >  a  we  have 


'a  +a    Va  +A    VA  +A 


Thus  at  no  point  ^   under  consideration,  (so  that  (C,  *  0) 
is  on  one  of  the  bounding  slits),  does  the  determinant 
vanish. 
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Now  note  that  Equation  (68)  says  that  the  point 
(x,y,t)  is  on  the  envelope  of  the  family  of  characteristic 
planes  given  by 

■KiK)    X  +  ti(C)y  =  t. 
(These  planes  are  characteristic,  for  we  have,  identically. 


P(X,ti)  =  l-(a^+A^)(X^+ti^)+a2AV(X^+ti^)  =  0. 

Thus  we  have  proved  assertion  (I4.),  that  the  boundary  curves 
in  the  plane  of  (x/t,  y/t)  represent  characteristic  surfaces. 
Assertion  (1)  follows  also;  for  the  line 

X(C)(Vt)  +  ti(C)(y7t)  =  1 

is  tangent  to  the  envelope  at 

(I)  =2^    (^)  =  2 
^t^    t  '  ^t^    t 

and  the  intercept  of  this  line  on   the  (x/t)  axis  is  at 
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(|>  =  xTTT  =  «  =  ^- 


To  move  assertion  (3),  we  expand  ti(z)  in  powers  of  z, 
We  have 


,  x.^^  ^  /  ^  /  v^-z  /a  ) 
z^i  =  i^l ^ ?,  2x  


J(l-z^/a^)(l- 
I     d-z^/c' 


where  in  this  formula  we  define  the  square  root  so  that  it 
Is  +1  at  z  =  Oo   Then  from 


X(x/t)  +  ix(y/t)  =  1 


we  obtain 


Expanding  in  powers  of  z,  we  have 

(f)  +  i(?)  [l+0(lzl^)]  =  z 
whence  assertion  (3)  follows. 
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Figure  I:     ENVELOPE   OF  HYDROMAGNTIC  PLANE  WAVE  FRONTS 


53  - 


Yt 


UNDISTURBED   fluid:     /),u,  B  =  0 


FRONT  OF  INCIDENT    WAVE 


HALF  PLANE  OBSTACLE 


DISTURBED     FLUID: 


Figure  2--     STATE  OF  FLUID   BEFORE    INCI  DENT  WAVE  STRIKES 
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NCIDENT  WAVE 


Figure  4: 


DIFFRACTION    OF  SLOW  SHOCK 
BY     HALF    PLANE 
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